A planetary system with gas giants and super-Earths around the nearby M dwarf GJ676A : Optimizing data analysis techniques for the detection of multi-planetary systems

G. Anglada-Escude and M. Tuomi, 'A planetary system via gas giants and super-Earths around the nearby M dwarf GJ 676A: Optimizing data analysis techniques for the detection of multi-planetary systems', Astronomy & Astrophysics, Vol. 548, A58 (2012). The version of record is available online at DOI: 10.1051/0004-6361/201219910 © ESO 2012 Published by EDP SciencesContext. Several M dwarfs are targets of systematical monitoring in searches for Doppler signals caused by low-mass exoplanet companions. As a result, an emerging population of high-multiplicity planetary systems around low-mass stars are being detected as well. Aims. We optimize classic data analysis methods and develop new ones to enhance the sensitivity towards lower amplitude planets in high-multiplicity systems. We apply these methods to the public HARPS observations of GJ≠676A, a nearby and relatively quiet M dwarf with one reported gas giant companion. Methods. We rederived Doppler measurements from public HARPS spectra using the recently developed template matching method (HARPS-TERRA software). We used refined versions of periodograms to assess the presence of additional low-mass companions. We also analysed the same dataset with Bayesian statistics tools and compared the performance of both approaches. Results. We confirm the already reported massive gas giant candidate and a long period trend in the Doppler measurements. In addition to that, we find very secure evidence in favour of two new candidates in close-in orbits and masses in the super-Earth mass regime. Also, the increased time-span of the observations allows the detection of curvature in the long-period trend. suggesting the presence of a massive outer companion whose nature is still unclear. Conclusions. Despite the increased sensitivity of our new periodogram tools, we find that Bayesian methods are significantly more sensitive and reliable in the early detection of candidate signals, but more work is needed to quantify their robustness against false positives. While hardware development is important in increasing the Doppler precision, development of data analysis techniques can help to reveal new results from existing data sets with significantly fewer resources. This new system holds the record of minimum-mass range (from Msini ~ 4.5 M⊕ to 5 Mjup) and period range (from P ~ 3.6 days to more than 10 years). Although all planet candidates are substantially more massive, it is the first exoplanetary system with a general architecture similar to our solar system. GJ 676A can be happily added to the family of high-multiplicity planetary systems around M dwarfs.Peer reviewe

Up to four planets around the M dwarf GJ 163 : Sensitivity of Bayesian planet detection criteria to prior choice

M. Tuomi and G. Anglada-Escude, 'Up to four planes around the M dwarf GJ 163: Sensitivity of Bayesian planet detection criteria to prior choice', Astronomy & Astrophysics, Vol. 556, A111 (2013), available online at DOI: 10.1051/0004-6361/201321174 © ESO 2013 Published by EDP SciencesContext. Exoplanet Doppler surveys are currently the most efficient means to detect low-mass companions to nearby stars. Among these stars, the light M dwarfs provide the highest sensitivity to detect low-mass exoplanet candidates. Evidence is accumulating that a substantial fraction of these low-mass planets are found in high-multiplicity planetary systems. GJ 163 is a nearby inactive M dwarf with abundant public observations obtained using the HARPS spectrograph. Aims. We obtain and analyse radial velocities from the HARPS public spectra of GJ 163 and investigate the presence of a planetary companions orbiting it. The number of planet candidates detected might depend on some prior assumptions. Since the impact of prior choice has not been investigated throughly previously, we study the effects of different prior densities on the detectability of planet candidates around GJ 163. Methods. We use Bayesian tools, i.e. posterior samplings and model comparisons, when analysing the GJ 163 velocities. We consider models accounting for the possible correlations of subsequent measurements. We also search for activity-related counterparts of the signals we observe and test the dynamical stability of the planetary systems corresponding to our solutions using direct numerical integrations of the orbits. Results. We find that there are at least three planet candidates orbiting GJ 163. The existence of a fourth planet is supported by the data but the evidence in favor of the corresponding model is not yet conclusive. The second innermost planet candidate in the system with an orbital period of 25.6 days and a minimum mass of 8.7 M· is inside the liquid-water habitable zone of the star. Conclusions. The architecture of GJ 163 system resembles a scaled-down Solar System in the sense that there are two low-mass planets on orbital periods of 8.7 and 25.6 days in the inner system, a possible slightly more massive companion on an intermediate orbit, and an outer sub-Saturnian companion at roughly 1 AU. The discovery of (yet) another planetary system with several low-mass companions around a nearby M-dwarf indicates that the high-multiplicity planetary systems found by the NASA Kepler mission around G and K dwarfs is also present (possibly even reinforced) around low-mass stars. © 2013 ESO.Peer reviewe

The origin of the excess transit absorption in the HD 189733 system: planet or star?

We have detected excess absorption in the emission cores of Ca II H&K during transits of HD 189733b for the first time. Using observations of three transits, we investigate the origin of the absorption, which is also seen in Hα and the Na I D lines. Applying differential spectrophotometry methods to the Ca II H and Ca II K lines combined, using respective passband widths of Δλ = 0.4 and 0.6 Å yields excess absorption of td = 0.0074 ± 0.0044 (1.7σ; Transit 1) and 0.0214 ± 0.0022 (9.8σ; Transit 2). Similarly, we detect excess Hα absorption in a passband of width Δλ = 0.7 Å, with td = 0.0084 ± 0.0016 (5.2σ) and 0.0121 ± 0.0012 (9.9σ). For both lines, Transit 2 is thus significantly deeper. Combining all three transits for the Na I D lines yields excess absorption of td = 0.0041 ± 0.0006 (6.5σ). By considering the time series observations of each line, we find that the excess apparent absorption is best recovered in the stellar reference frame. These findings lead us to postulate that the main contribution to the excess transit absorption in the differential light curves arises because the normalizing continuum bands form in the photosphere, whereas the line cores contain a chromospheric component. We cannot rule out that part of the excess absorption signature arises from the planetary atmosphere, but we present evidence which casts doubt on recent claims to have detected wind motions in the planet's atmosphere in these data

Color Difference Makes a Difference: Four Planet Candidates around τ Ceti

The removal of noise typically correlated in time and wavelength is one of the main challenges for using the radial-velocity (RV) method to detect Earth analogues. We analyze τ Ceti RV data and find robust evidence for wavelength-dependent noise. We find that this noise can be modeled by a combination of moving average models and the so-called "differential radial velocities." We apply this noise model to various RV data sets for τ Ceti, and find four periodic signals at 20.0, 49.3, 160, and 642 days, which we interpret as planets. We identify two new signals with orbital periods of 20.0 and 49.3 days while the other two previously suspected signals around 160 and 600 days are quantified to a higher precision. The 20.0 days candidate is independently detected in Keck data. All planets detected in this work have minimum masses less than 4M⊕ with the two long-period ones located around the inner and outer edges of the habitable zone, respectively. We find that the instrumental noise gives rise to a precision limit of the High Accuracy Radial Velocity Planet Searcher (HARPS) around 0.2 m s−1. We also find correlation between the HARPS data and the central moments of the spectral line profile at around 0.5 m s−1 level, although these central moments may contain both noise and signals. The signals detected in this work have semi-amplitudes as low as 0.3 m s−1, demonstrating the ability of the RV technique to detect relatively weak signals

AD Leonis: Radial Velocity Signal of Stellar Rotation or Spin–Orbit Resonance?

AD Leonis is a nearby magnetically active M dwarf. We find Doppler variability with a period of 2.23 days, as well as photometric signals: (1) a short-period signal, which is similar to the radial velocity signal, albeit with considerable variability; and (2) a long-term activity cycle of 4070 ± 120 days. We examine the short-term photometric signal in the available All-Sky Automated Survey and Microvariability and Oscillations of STars (MOST) photometry and find that the signal is not consistently present and varies considerably as a function of time. This signal undergoes a phase change of roughly 0.8 rad when considering the first and second halves of the MOST data set, which are separated in median time by 3.38 days. In contrast, the Doppler signal is stable in the combined High-Accuracy Radial velocity Planet Searcher and High Resolution Echelle Spectrometer radial velocities for over 4700 days and does not appear to vary in time in amplitude, phase, period, or as a function of extracted wavelength. We consider a variety of starspot scenarios and find it challenging to simultaneously explain the rapidly varying photometric signal and the stable radial velocity signal as being caused by starspots corotating on the stellar surface. This suggests that the origin of the Doppler periodicity might be the gravitational tug of a planet orbiting the star in spin–orbit resonance. For such a scenario and no spin–orbit misalignment, the measured v sin i indicates an inclination angle of 15°̣5 ± 2°̣5 and a planetary companion mass of 0.237 ± 0.047 M Jup

High-cadence spectroscopy of M-dwarfs – II. Searching for stellar pulsations with HARPS

Stellar oscillations appear all across the Hertzsprung–Russell diagram. Recent theoretical studies support their existence also in the atmosphere of M dwarfs. These studies predict for them short periodicities ranging from 20 min to 3 h. Our Cool Tiny Beats (CTB) programme aims at finding these oscillations for the very first time. With this goal, CTB explores the short time domain of M dwarfs using radial velocity data from the High Accuracy Radial velocity Planet Searcher (HARPS)-European Southern Observatory and HARPS-N high-precision spectrographs. Here we present the results for the two most long-term stable targets observed to date with CTB, GJ 588 and GJ 699 (i.e. Barnard's star). In the first part of this work we detail the correction of several instrumental effects. These corrections are especially relevant when searching for subnight signals. Results show no significant signals in the range where M dwarfs pulsations were predicted. However, we estimate that stellar pulsations with amplitudes larger than ∼0.5 m s−1 can be detected with a 90 per cent completeness with our observations. This result, along with the excess of power regions detected in the periodograms, opens the possibility of non-resolved very low amplitude pulsation signals. Next generation more precise instrumentation would be required to detect such oscillations. However, the possibility of detecting pulsating M-dwarf stars with larger amplitudes is feasible due to the short size of the analysed sample. This motivates the need for completeness of the CTB survey

Analysis of apsidal motion in eclipsing binaries using TESS data: I. A test of gravitational theories

The change in the argument of periastron of eclipsing binaries, i.e., the apsidal motion caused by classical and relativistic effects, can be measured from variations in the difference between the time of minimum light of the primary and secondary eclipses. Poor apsidal motion rate determinations and large uncertainties in the classical term have hampered previous attempts to determine the general relativistic term with sufficient precision to test General Relativity predictions. As a product of the TESS mission, thousands of high-precision light curves from eclipsing binaries are now available. Using a selection of suitable well-studied eccentric eclipsing binary systems, we aim to determine their apsidal motion rates and place constraints on key gravitational parameters. We compute the time of minimum light from the TESS light curves of 15 eclipsing binaries with precise absolute parameters and with an expected general relativistic contribution to the total apsidal motion rate greater than 60%. We use the changing primary and secondary eclipse timing differences over time to compute the apsidal motion rate, when possible, or the difference between the linear periods as computed from primary and secondary eclipses. For a greater time baseline we carefully combine the high-precision TESS timings with archival reliable timings. We determine the apsidal motion rate of 9 eclipsing binaries, 5 of which are reported for the first time. From these, we are able to measure the general relativistic apsidal motion rate of 6 systems with sufficient precision to test General Relativity for the first time using this method. This test explores a regime of gravitational forces and potentials that had not been probed earlier. We find perfect agreement with the theoretical predictions, and we are able to set stringent constraints on two parameters of the parametrised post-Newtonian formalism.Comment: Accepted for publication in A&

High-cadence spectroscopy of M dwarfs – I. Analysis of systematic effects in HARPS-N line profile measurements on the bright binary GJ 725A+B

Understanding the sources of instrumental systematic noise is a must to improve the design of future spectrographs. In this study, we alternated observations of the well-suited pair of M-stars GJ 725A+B to delve into the sub-night High Accuracy Radial Velocity Planet Searcher for the Northern hemisphere (HARPS-N) response. Besides the possible presence of a low-mass planet orbiting GJ 725B, our observations reveal changes in the spectral energy distribution (SED) correlated with measurements of the width of the instrumental line profile and, to a lower degree, with the Doppler measurements. To study the origin of these effects, we searched for correlations among several quantities defined and measured on the spectra and on the acquisition images. We find that the changes in apparent SED are very likely related to flux losses at the fibre input. Further tests indicate that such flux losses do not seriously affect the shape of the instrumental point spread function of HARPS-N, but identify an inefficient fitting of the continuum as the most likely source of the systematic variability observed in the full width at half-maximum. This index, accounting for the HARPS-N cross-correlation profiles width, is often used to decorrelate Doppler time series. We show that the Doppler measurement obtained by a parametric least-squares fitting of the spectrum accounting for continuum variability is insensitive to changes in the slope of the SED, suggesting that forward modelling techniques to measure moments of the line profile are the optimal way to achieve higher accuracy. Remaining residual variability at ∼1 m s−1 suggests that for M-stars Doppler surveys the current noise floor still has an instrumental origin

Correcting for chromatic stellar activity effects in transits with multiband photometric monitoring: Application to WASP-52

The properties of inhomogeneities on the surface of active stars (i.e. dark spots and bright faculae) significantly influence the determination of the parameters of an exoplanet. The chromatic effect they have on transmission spectroscopy could affect the analysis of data from future space missions such as JWST and Ariel. To quantify and mitigate the effects of those surface phenomena, we developed a modelling approach to derive the surface distribution and properties of active regions by modelling simultaneous multi-wavelength time-series observables. By using the StarSim code, now featuring the capability to solve the inverse problem, we analysed $\sim$ 600 days of BVRI multiband photometry from TJO and STELLA telescopes exoplanet host star WASP-52. From the results, we simulated the chromatic contribution of surface phenomena on the observables of its transits. We are able to determine the relevant activity parameters of WASP-52 and reconstruct the time-evolving longitudinal map of active regions. The star shows a heterogeneous surface composed of dark spots with a mean temperature contrast of $575\pm150$ K with filling factors ranging from 3 to 14 %. We studied the chromatic effects on the depths of transits obtained at different epochs with different stellar spot distributions. For WASP-52, with peak-to-peak photometric variations of $\sim$7 % in the visible, the residual effects of dark spots on the measured transit depth, after applying the calculated corrections, are about $10^{-4}$ at 550 nm and $3\times10^{-5}$ at 6$\mu$m. We demonstrate that by using contemporaneous ground-based multiband photometry of an active star, it is possible to reconstruct the parameters and distribution of active regions over time, and thus, quantify the chromatic effects on the planetary radii measured with transit spectroscopy and mitigate them by about an order of magnitude.Comment: 17 pages, 12 figure